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Data from e ! e annihilation into hadrons collected by the JADE experiment at centre-of-mass 
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a s (M z o) = 0.1286 ± 0.0007(stat.)± 0.0011 (exp.) ± 0.0022(had.) ± 0.0068 (theo.), 
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a s (M z o) = 0.1169 ± 0.0004(stat.)± 0.0012(exp.) ± 0.0021 (had.) ± 0.0007(theo.), 
both in agreement with the world average. 
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1. Introduction 

The annihilation of an electron and a positron into hadrons allows precise tests of Quantum 
Chromodynamics (QCD). We measured the first five moments of event shape observables and the 
four-jet rate in e" e annihilation and compared the data to predictions by Monte Carlo Models and 
by perturbative QCD. From the comparison of the data with the theory we extracted the strong cou¬ 
pling as- A more detailed description of this analysis can be found in |[I]]. We used data collected 
by the JADE experiment [j2|] in the years 1979 to 1986 at the PETRA e 1 e collider at DESY at 
six center-of-mass energies of = 14.0,22.0,34.6,35.0,38.3 and 43.8 GeV. The total integrated 
luminosity was 195 pb 1 corresponding to about 40000 events. Large samples of Monte Carlo sim¬ 
ulated events were used to correct the data for experimental acceptance, resolution, backgrounds 
and fragmentation effects. The process e + e " —► hadrons was simulated using PYTHIA and cor¬ 
responding samples using F1ERWIG and ARIADNE were used for systematic checks. The event 
selection for this analysis aimed to identify hadronic event candidates and rejected events with a 
large amount of energy emitted by initial state radiation (ISR), leptonic final states and two photon 
events. Events from the process e 1 e ~ —> bb were considered as background, since especially at low 
center-of-mass energy the large mass of the b quarks and of the subsequently produced B hadrons 
will influence the measurement. Therefore the contribution from expected bb events was subtracted 
[3], The effects of detector acceptance and resolution and of residual ISR were accounted for by 
a multiplicative correction procedure. In order to confront the theoretical prediction valid for par- 
tons with the data the theoretical prediction was corrected for hadronization with PYTHIA. Several 
sources of possible systematic uncertainties were studied. The experimental uncertainties evaluate 
deficiencies in the reconstruction of the data and the modelling of the detector. The uncertainties 
associated with the correction for hadronization effects were assessed by using alternative models 
(HERWIG and ARIADNE instead of the default PYTHIA). The theoretical uncertainty associated 
with missing higher order terms in the theoretical prediction, was assessed by setting the renormal¬ 
ization scale Xu to 0.5 and 2 instead of the default value one. The results from the fits at the various 
center-of-mass energies were combined in order to determine a single value of cxs(Myn). 

2. Event Shape moments 

The properties of hadronic events may be characterized by event shape observables. The fol¬ 
lowing event shapes were considered in this analysis: Thrust, C-Parameter, Heavy Jet Mass, Jet 
Broadening Observables Bj and B\\ and the transition value between 2 and 3 jets >’23 using the 
Durham scheme. The 72 th, n = 1,2,... ,5 moment of the distribution of an event shape observable 
y was defined by (y n ) = f^ mwc y n ^^dy , where y max was the kinematically allowed upper limit of 
the observable. The measurement involved a full integration over the available phase space. Com¬ 
parisons of QCD predictions were thus complementary to tests of the theory using the differential 
distributions. In the case of the moments of event shape observables the QCD predictions were 
obtained by numerical integration of the QCD matrix elements using the program EVENT2 |)4]. 
A x 2 value was minimized with respect to a,s(M z o) for each moment (/') separately. The fit re¬ 
sults are shown in figure [I]. The fit (Mu) did not converge and therefore no result is shown. We 
observed values of ^ 2 /d.o.f. of G( I); the fitted QCD predictions including the running of as were 
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Figure 1: The figure on the left shows the measurements of as(M z o) using fits to moments of six event 
shape observables. The inner error bars represent statistical errors, the middle error bars include experimental 
errors and the outer error bars show the total errors. The dotted line indicates the weighted average described 
in the text; only the measurements indicated by solid symbols were used for this purpose. The figure on the 
right shows the values for as using the four-jet rate at the various energy points. The errors show the 
statistical (inner part) and the total error. The full and dash-dotted lines indicate the current world average 
value of as(M z o) |R|. The results at \/s = 34.6 and 35 GeV have been combined for clarity. The results 
from ALEPH [g] and OPAL [01 (preliminary) are shown as well. 



thus consistent with our data. However, we fond that for ((1 — T) n ), (C n ) and (B'j) the fitted values 
of a,s(M z o) increased steeply with the order n of the moment used. This effect was in clear cor¬ 
relation with increasing ratio of the NLO and LO coefficient with moment n for ((1 — T) n ), (C n ") 
and We considered only those results for which the NLO term was less than half the LO 

term, namely (1 — T), (C), ( Bj), (B and ((V23)''), n = 1,... ,5 and (M^), n = 2,... ,5; i.e. results 
from 17 observables in total. The purpose of this requirement was to select observables with an 
apparently converging perturbative prediction. The statistical correlations between the 17 results 
were determined using Monte Carlo simulation at the hadron level. The result of the combination 
was 


a s (M z 0 ) = 0.1286± 0.0007(stat.) ±0.0011(exp.) ±0.0022(had.) ± 0.0068 (theo.) , 

above but still consistent with the world average value of as(M z o ) = 0.1182 ±0.0027 [6]. Com¬ 
bining only the fit results from (1 — T), (C ), (Bj), (Bw)> {>’23 ) and (Mg) yields a value of 

as(M z o) = 0.1239±0.0001 (stat.) ±0.0008(exp.) ±0.0009(had.) ±0.0059(theo.) . 

The slightly smaller error of as reflects the fact that the lower order moments were less sensitive to 
the multi-jet region of the event shape distributions. This leads to a smaller systematic uncertainty. 
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2.1 Four-Jet Rate 

Besides event shape observables jet rates reflect the parton structure of the event. The four- 
jet rate R 4 (y cu t) is the fraction of four-jet events as a function of v cu t- the jet-resolution parameter 
determined in the Durham scheme. In QCD the fraction of four-jet events R 4 is predicted in NLO as 
a function of the strong coupling as- The prediction was calculated by the program DEBRECEN 
2.0 [§] combined with an all order resummation using the “modified R-matching” scheme. In 
the case of jet rates a single event usually contributed to several y cut points in the four-jet rate 
distribution and for this reason the data points were correlated. The complete covariance matrix was 
determined and used in the y 2 fit for the extraction of as- The y 2 value was minimized with respect 
to as for each center-of-mass energy separately. The fit ranges covered the decreasing parts of the 
distributions at large y cu t, where the perturbative QCD predictions were able to adequately describe 
the data corrected for hadronization. We found that the fit result from the 14 GeV data had large 
hadronization and experimental uncertainties because the corresponding corrections were large and 
not well known at this energy. We therefore chose not include this result in the combination. The 
result of the combination using all results with \fs > 22 GeV was 

a s (M z o) =0.1169± 0.0004(stat.) ±0.0012(exp.) ±0.0021(had.) ±0.0007(theo.) , 

consistent with the world average value of as(M z o) = 0.1182 ± 0.0027 [^]. The results at each 
energy point are shown in figure [I] and compared with the predicted running of as based on the 
world average value. 

3. Summary 

In this note we present preliminary measurements of the four-jet rate and the first five moments 
of event shape observables at center-of-mass energies between 14 and 44 GeV using data of the 
JADE experiment. The value of the strong coupling is determined to be as(M z o) = 0.1169 ± 
0.0026(tot.) using the four-jet rate and as(M z o) = 0.1286±0.0072(tot.) using moments of event 
shape observables. 
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